Elusive gauge-invariant fermion propagator in QED-like effective theories: round II 
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We comment on the recent attempt by M. Franz et al [Q| to further justify their earlier calculation 
of the gauge-invariant electron propagator in the context of the QED3 theory of the pseudogap 
phase in cuprates First, we use the method of "reductio ad absurdum" to demonstrate the 
inconsistency of the argument offered in jl| and then present a direct calculation of the disputed 
fermion amplitudes, thus unequivocally proving that the previously proposed form of the electron 
propagator: 1) does exhibit a negative anomalous dimension, as pointed out in 2) is different 
from the so-called Brown's function whose anomalous dimension turns out to be positive when 
computed in a covariant gauge. Lastly, we conjecture that the true physical electron propagator 
(whose exact form still remains to be found) may feature a " super- Luttinger" behavior characterized 
by a faster than a power-law decay Gphys{x) oc exp(— const ln^(A[x|)). 



In their recent note Q], M. Franz et al continued to advocate the naive "stringy ansatz" with the inserted Wilson 
line-like factor taken along the straight path T between the end points x and y 

Goix - y) 0|V'(x) exp(-i ^ A^(z)dz^)^(y)|0 > (1) 

as a viable candidate for the gauge invariant propagator of physical electrons in the QED3 theory of the pseudogap 
phase in cuprates In this theory, the electron operators are related by virtue of a singular gauge transformation 
^f(x) = e'^(^^°°f^V(a;) to the A^-flavored La grangian fermion variables (in the theory of the pseudogap phase of Ref. 
= 2) which are governed by the QED3 action 

_ f ^ _ 1 

S[^, ^P,A]^ / dz[^ i'fiijA, + %A,,)^jf + -^id,A, - d,A^)^] (2) 

f=i 9 

In the physically interesting strong coupling regime of momenta q ^ A = Ng^ the effect of fermion polarization is 
described by the renormalized gauge field propagator 



D^,{q) =< Q\A^{q)A,{-q)\Q >= ^V^IS^"" + (A - 1)^] (3) 

8 (7^ 



which is parametrized by A in the class of covariant (generalized Lorentz) gauges. 

In particular, the authors of [|j claimed that the amplitude (1) exhibits a physically sensible behavior characterized 
by a positive anomalous dimension. This conclusion was drawn solely on the basis of a seemingly convincing argument 
that the amplitude (1) can be identified with another function known as the Brown's construct jj] 

G,(x-y) = < m-my)\o > .4. 

whose anomalous dimension can be deduced rather straightforwardly from the ratio between the wave function renor- 
malization factor determining the anomalous dimension of the ordinary (gauge variant) fermion propagator 

_ < Om^my)\0 > ^ ^|)(4/3.^.)(.-3.) (5) 

G{x - y) 

(hereafter G{x — y) stands for the bare propagator) and the Gaussian average of the Wilson line 

<eM^ j^A^dz^) >=exphi^dz5'^dz2^7^^,(zi-Z2)](x (A|a;-j/|)(4/-'^)(2-A), (6) 

thus resulting in the overall positive anomalous dimension |^ 
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which is free of the gauge parameter A, thus creating the impression that Eq.(4) represents a truly gauge invariant 
function. If this were the case, then one would indeed be able to identify Eqs.(l) and (4) by simply noticing that the 
two functions appear to coincide in the axial (Fock-Shwinger) gauge (a; — y)^A^(z) = where the Wilson line factor 
equals unity. 

In their recent follow-up note [Q, M. Franz ct al attempted to further support these conclusions by invoking a 
formal textbook argument |^ which allows one to switch between different gauges when computing a gauge invariant 
quantity such as, e.g., partition sum. In fact, it is a confidence in the (already established) gauge invariance of an 
amplitude in question that makes this argument a meaningful statement, while, in the absence of such, this whole 
argument becomes largely irrelevant. 

In our earlier note we argued that the argument of Ref. does have a caveat, for the Brown's function defined 
as a ratio of the two different amplitudes can not be truly gauge invariant, for the gauge fields in the numerator and 
denominator of Eq.(4) may transform totally independently of one another (A^'^ Aj^'^ + d^f^'^), thus resulting in 

the overall phase factor e*'^'^'^)"*-''^'^^)"*-''^^^^^*'^^'^^) which only vanishes for f^{x) — f^{x). 

Albeit not exactly gauge invariant, the Brown's function is, nevertheless, independent of the gauge parameter A in 
the particular class of the covariant gauges, as manifested by Eq.(7). What invalidates the attempt to identify Eqs.(l) 
and (4) despite their coinciding with each other in the axial gauge is the fact that the latter gauge is not a covariant 
one. Therefore, the surrogate Eq.(4) can not substitute for the amplitude (1) whose dimension was previously found 
to be negative 

thus making Eq.(l) a rather poor candidate for the physical electron propagator. 

In fact, if the conjectured identity between Eqs.(l) and (4) proved to be true, it would also hold for any function 
defined as 

< 0|V>(x)exp(z(e- l)/r A^(z)rfz,.)|0 > ?(y)|0 > 
^^"^ ^' <0|exp(z$/j,A^(z)dz^)|0> ' 

because: 1) Eq.(9) appears to be seemingly gauge invariant to the same extent as Eq.(4) (that is, provided that one 
uses the same gauge transformation in both numerator and denominator); 2) it also coincides with both Eqs.(l) and 
(4) when computed in the axial gauge (observe that Eqs.(l) and (4) correspond to ^ = and ^ = 1, respectively). 
The parameter ^ should not be confused with the gauge parameter A (for a reader's convenience, we point out that 
in Refs. |^,^ the gauge parameter is denoted as ^, while our parameter ^ has no counterpart). 

It can be readily seen, however, that the alleged ^-independence of Eq.(9) would, in fact, be too much to wish for. 
To this end, we represent Eq.(9) as a functional average (hereafter denoted by brackets <...>) 

<gix,y\A)e^(i-m^-y\^)> 

Gd^-y) = < e^m.,y\A) > ' (10) 

where 9{x,y\A) = Jp A^(z)dz^ and Q{x,y\A) — l/{id + A) stands for the inverse Dirac operator, over different gauge 
field configurations with the weight determined by the effective action 

Seff[A] = ljd^j dyA^ix)D;^ix - y)A,iy) (11) 

where, in accord with all the previous work on the subject, we neglect higher order (non-Gaussian) corrections 
produced by fermion polarization, thereby focusing on the leading terms in the 1/N expansion. 

It can be readily seen that, if Eq.(9) were indeed independent of ^, one would be able to establish a number of 
clearly improbable identities, including 

< e(a;,y|A)e''(«-i)^("'^l^) >=< ^(x, A)e-*(«+i)*("'^l^) > (WRONG) (12) 

which stems from equating (9) to its value obtained for — ^ and noticing that in a covariant gauge (and in the leading 
1/N approximation) one has < e'?^^^'?''^) >=< e-^«^(^'2'l^) > (see (6)). 

In particular, if valid, Eq.(12) would have implied that the ordinary (gauge- variant) fermion propagator < 
0\ip{x)ilj{y)\0 >=< Qix, y\A) > coincides with the amplitude < t/(x, i/|A)e~^'^^^'^l^'' > not only in the axial (where it 
does hold) but also in an arbitrary, including any covariant, gauge. 
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And it gets even better: by differentiating Eq.(9) with respect to ^ and putting £, = one finds that the alleged 
^-independence (hence, the requirement dG^{x — y)ld^\^=Q = 0) imposes the condition 

I < g{x,y\A)e{x,y\A)e-''''-^-y\'^^ >= {WRONG) (13) 

The validity of (13) can be ascertained by simply computing this amplitude, for which purpose it would be most 
convenient to use the exact quantum mechanical (i.e., single-particle) path integral representation of the inverse Dirac 
operator G{x, y\A) (also, see |l0| for an asymptotically exact non-perturbative calculation of the gauge-invariant 

fermion propagator in the case of a static gauge field which one encounters in such problems as the effect of vortex 
disorder on quasiparticle properties of c?-wave superconductors or that of dislocations in layered graphite). However, 
here we choose a more traditional approach, for we recognize that the more powerful method devised in |^J^ and 
further advanced in PJTo[| may still be somewhat less familiar. 

In fact, it suffices to compute (13) to first order in 1/A^ by expanding both the exponential factor e"'^'^^^^''^-' = 
1 — i9{x,y\A) + . .. and the inverse Dirac operator l/{id + A) ^ 1/id — {1 / id) A{1 / id) + . . . in powers of Af^{z) 
(apparently, no cancellation can possibly occur between different 1/A^ order terms) 

I < g{x, y\A)0{x, y|A)e-'^("'^l^) >= G{x - y) dz^ ^ dz^D^,{z, - z^)^ 

f 16 

dziG{x - zi)j^G{zi - y) / dz!^D,,,{zi - z^) j-G(a; - y) ln(A|x - y|) ^ {Q.E.D.) (14) 

Jr ""A 

where we used the following Eqs.(17) and (18) (see below) after scaling out the factors (1 — 2^) and (1 — ^), respectively. 
Taken at its face value, the result (14) implies that Gq(x — y) ^ Gi{x — y), in complete accord with the proverbial 
moral: " the average of a product is not necessarily equal to a product of the averages" . 

In fact, the sought-after identity between the different functions G^{x) can approximately hold only in the case 
of massive fermions and only in the vicinity of the mass shell (|p^ — <C m^), resulting in their common (rather 
trivial) long-distance behavior G^{x) oc e"™'^'. In this regime, the leading functional dependence on the gauge field 
reduces to the eikonal phase factor g{x,y\A) ^ g^^i^^vl^) and it is precisely this property that the Brown's 

function was designed to make use of in the first place Q . 

However, in the massless case such a regime is absent altogether, and the long-distance behavior is solely controlled 
by the ultra-violet (UV) anomalous dimensions (see for a more extensive discussion of this potentially confusing 
issue). 

Having completed our proof, we now present a direct calculation of (x) for an arbitrary ^ in a generic covariant 
A-gauge. After expanding (9) to the second order in Ap(z), we find three different kinds of correction terms 

G^{x -y)~ G{x ~y)^^ j dzidz2 < G{x - zi)A{zi)G{zi - z2)A{z2)G{z2 - y) > + 

+ ^-^^^G{x -y)< (Oix, y\A))' > +(1 - ^ J dz < G{x - z)A{z)G{z ~ y)e{x, y\A) > (15) 

First of these terms corresponds to the lowest order self-energy correction to the ordinary Green function 

SiG^ix -y) = -J dzidz2Gix ~ zi)-i^D^,{zi - 22)7.^(^2 - y) - G{x - Z/)^^(2 - 3A) ln(A|x - y\) (16) 

thus reproducing the anomalous dimension associated with the wave function renormalization factor (5). 

The second type of corrections originates from the expansion of the Wilson lines inserted to the numerator and 
denominator of Eq.(9) 

52G^{x~y) = ^—^^^^G{x-y) [ dz^ [ dz^D,,{z^ - Z2) = Gix - y)-^il ~ 2^(2 - X)HA\x ~ y\) (17) 



r Jr 



which, too, can be readily exponentiated and yields the anomalous dimension equal to that given by Eq.(6) with the 
extra factor (1 — 2^). 

Lastly, there is a mixed term which stems from the first order expansion of both G{x,y\A) and the Wilson line. 
Unlike Eqs.(16) and (17), this contribution has not been discussed in the previous studies, and, therefore, we spell its 
evaluation out in all the details 

hG^ix - y) = (1 - f) J dziG{x - zi)7pG(zi ~ v) J dz2D^^{zi - Z2) ~ 
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(1 - y d7^i[G{x - yh^,G{zi -y) + G{x - zi)j^G{x - y)] dzl^D^^zi - Z2) 



(1 - OG(x - y)7;, / dz^ / dzi ^[A<5^, + 2(l-A) 



n^N' ^'^'-^ ""rlr \l |^l-^2P^ ^ ' \ZI-Z2\^ 



= G(x-y)^(l-0Aln(A|x-2/|) (18) 



Eq.(18) was obtained with logarithmic accuracy, and in the course of this calculation we used the 3D real space 

propagators 



and the following D-dimensional integrals 



/ 



and 



/ 



^ 2t,d/2 2y^5P^ -yH^^ -y^d'^P y'^y^y^ 

3(Z)-l)r(i5/2)^ |y|c-i ^^-^^ 



Combining Eqs.(16), (17), and (18) together we finally obtain the anomalous dimension of G^{x) in the form 

# = 3^(3«-2) (22) 

Thus, despite the fact that the gauge parameter A cancels out, as expected, the functions G^{x) are starkly different 
for different ^. In particular, Eqs.(7) and (8) follow from (22) for ^ = 1 and ^ = 0, respectively. 

Anticipating possible objections that the above results might have been different, should we have resorted to the 
dimensional rcgularization of our divergent integrals, we have also confirmed that in D ^ 3 dimensions Eqs.(16), 
(17), and (18) are all proportional to dzf'iz - /\z - y\^-^ = {A\x - y\f-^ / (3 - D) ln(A|x-y|). 

For the sake of completeness, here wc also present the results for the weak coupling QED4 which demonstrate that 
the situation in 3D is not at all exceptional. Instead of Eqs.(16),(17), and (18) we now get 

5^G^{x -y) = -Gix - y)-£^Xln{A\x - y\) (23) 

52G^ix -y)= G{x - y)-^(l - 2e)(3 - A) ln(A|x - y\) (24) 
6sG^{x -y) = G{x - y)-^{l - 0\\n{K\x - y\) (25) 



Instead of Eq.(19) we use 



and because of a different power-counting the 4D calculation involves the integrals 

/dx ^- = ^ 

J \x\D\x-y\D-^ T{L 



,D/2 ya 

and 



(27) 
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'''' mx-v\^ - 2{D~2)T{D/2Y jy^-^ + l^;^ ^ 

Combining Eqs. (23), (24), and (25) together we finally obtain the total anomalous dimension of G^{x) in the 3D 
case 

^f = g^(2^-l) (29) 

Interestingly enough, the values of Eq.(29) obtained for ^ = and ^ = 1 differ only in their sign. 

We emphasize that the positive values of i]^^''^^ pertaining to the (not exactly gauge-invariant) function 6*1(0;) 
computed in the covariant gauge does not fix the problem with the negative anomalous dimension rj^^''^^ of the 
conjectured form of the physical electron propagator which is given by the (exactly gauge-invariant) function Gq{x). 
In light of these findings, in Refs. we suggested to actively explore alternate proposals for the physical electron 
propagator which, due to the nature of the relationship between the electrons and the Lagrangian fermions imple- 
mented through a singular gauge transformation, can only be given by a single gauge field average, not a ratio of 
such. 

To this end, in jsplwe discussed the "dressed charge" propagator constructed in the context of the conventional 
(i.e., massive) Q£^Z34 jlll] . The authors of Ref. showed that the Fourier transform of this gauge invariant amplitude 



- ,) =< oi^(.).xp| /;^:|;;;;-;;;:^;;-^;;> ,^.w.(.)io > (30) 



where m'^ = (1, 0) and — (0, v) features a simple pole-like behavior at the single point p'* = "^{1, v)/y/l — on the 
mass shell which corresponds to a charge moving with a velocity v. The UV anomalous dimension of Gr,{x) computed 
inQ 

^.'^=-f^[3 + 2-lni-^] (31) 

OTT^ V 1 + f 

remains positive for any v and increases logarithmically towards infinitely high values as v ^ 1. 

Albeit not being immediately applicable to the case m = 0, the calculation carried out in Ref. suggests that, 
as one proceeds beyond the first order, the logarithmic growth of (31) gets cut off at max{\v — 1|, 1/A|x|). Therefore, 
it is not totally inconceivable that the massless counterpart of Eq. (30) may exhibit a faster than a power-law decay 

Gphys{x) cx exp(— const ln^(A|a;|)), (32) 

where, depending on the dimension, the constant is proportional to either 1/N or g^ , thus placing the effective QED- 
like theories of condensed matter systems into the class of " super-Luttinger" models, alongside the ID metals with 
unscreened Coulombic interactions where Gphys{x) oc exp(— consi ln'^^^(A|2:|)). 

In light of such a possibility, the previous attempts to discover the Luttinger-like behavior in the pseudogap phase 
of the cuprates may need to be prepared to handle a potentially much stronger suppression of the physical 
amplitudes in order to reconcile the predictions of the QED3 theory with the available photoemission, tunneling, and 
other data. 

We conclude by stressing that, arguably, the problem of constructing the gauge invariant fermion propagator in 
massless QED, thus far, has received a lesser attention than it deserves. The intricacy of the related calculations 
indicates that this problem really needs to be settled before one can start drawing solid (instead of wishful) conclusions 
about the true behavior in the QED-like as well as other gauge field models, including non-abelian and discrete 
symmetry (say, Zn) ones. A successful completion of this task is likely to require some new ideas and/or potentially 
cumbersome calculations. Quoting from yet another textbook, in addition to the already mentioned ones jj,^, simply 
wouldn't do. 

The author acknowledges valuable email communications with V. Gusynin and I. Herbut. This research was 
supported by the NSF under Grant No. DMR-0071362. 
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